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Abstract 

y—i The advent of transformation optics and metamaterials has made possible devices producing extreme 

. !_H effects on wave propagation. Here we give theoretical designs for devices, Schrodinger hats, acting as 

X 

;_( invisible concentrators of waves. These exist for any wave phenomenon modeled by either the Helmholtz 



or Schrodinger equations, e.g., polarized waves in EM, pressure waves in acoustics and matter waves in QM, 
and occupy one part of a parameter space continuum of wave-manipulating structures which also contains 
standard transformation optics based cloaks, resonant cloaks and cloaked sensors. For EM and acoustic 
Schrodinger hats, the resulting centralized wave is a localized excitation. In QM, the result is a new charged 
quasiparticle, a quasmon, which causes conditional probabihstic illusions. We discuss possible solid state 
implementations . 

PACS numbers: 



Transformation optics and metamaterials have made possible devices producing effects on wave 
propagation not seen in nature, including invisibility cloaks for electrostatics [[Dl2l, electromag- 
netism (EM) flSHll, acoustics [I6H1 and quantum mechanics (QM) OH; field rotators HOll; EM 
wormholes [[TT]|: and illusion optics [[121, among many others. The purpose of this paper is to give 
theoretical designs for devices, which we refer to as Schrodinger hats, acting as invisible concen- 
trators, reservoirs and amplifiers for waves. Schrodinger hats exist for any wave phenomenon 
modeled by either the Helmholtz or Schrodinger equation, whether in EM, acoustics or QM. 
Schrodinger hats (SH) occupy one part of a parameter space continuum of wave-concentrating 
structures which also contains standard transformation optics based cloaks and cloaked sensors 
[|T3ll . For EM and acoustic SH, the resulting centralized wave is a localized excitation, which may 
be super-wavelength in scale; in QM, the SH produces a new quasiparticle, a quasmon. Acoustic 
and EM Schrodinger hats require negative index materials, while highly oscillatory potentials are 
needed for QM hats. A SH seizes a large fraction of an incident wave, holding and amplifying it 
as a quasmon, while contributing only a negligible amount to scattering. These devices are con- 
sistent with the uncertainty principle, and we illustrate the concept by a QM version of three card 
monte. While a quantum Schrodinger hat is invisible to one -particle scattering, we show by effec- 
tive potential theory f\^ that a SH acts as an amplifier of two-particle Coulomb interactions. Such 
amplifiers may be useful for quantum measurement and information processing. The similar yet 
less demanding acoustic and EM hats offer equivalent effects, but existing metamaterials [1151 - 1181 
make these designs more immediately realizable, allowing verification and further exploration of 
Schrodinger hats and quasmons. 

There are a number of 'paradoxes' in which the laws of quantum mechanics imply results that 
conflict with our intuition [[T9l . In this spirit, here we show that the behavior of matter waves, 
as governed by Schrodinger's equation, combined with the virtual space/physical space paradigm 
of transformation optics, allows one to manipulate conditional probabilities in QM and create 
quantum illusions, in which observed locations of particles differ from their actual values. 

Ideal (perfect) 3D quantum invisibility cloaks at fixed energy E are based on the behavior 
of solutions to Schrodinger equations, with specific potentials and singular, inhomogeneous and 
anisotropic mass density [9[. These are mathematically equivalent, via a Liouville gauge transfor- 
mation, to Helmholtz equations which also allow for cloaking in scalar optics [[3l[20|[ and acoustics 
[|71 [HI |2T1. Realizing a QM cloak would be challenging, due to the extreme material parame- 
ters required [9[. We have previously described approximate QM cloaks, avoiding extreme and 




FIG. 1: A quasmon inside a Schrodinger hat. The real part of the effective wave function t/j^^Jx, y, z) 
at the plane z = when a plane wave is incident to a SH potential. By varying the design parameters, the 
concentration of the wave inside the cloaked region can be made arbitrarily strong and the scattered field 
arbitrarily small. The matter wave is spatially localized, but conforms to the uncertainty principle, with the 
large gradient of ip, visible as the steep slope of the central peak, concentrating the momentum in a spherical 
shell in p-space. 



anisotropic parameters but nevertheless acting with arbitrary cloaking effectiveness Il22ll23l . If a 
matter wave is incident to such a potential, the scattered wave can be made as small as desired. 
Analysis of approximate QM cloaks revealed a difficulty: the wave vanishes inside the cloak un- 
less the cloak supports an almost trapped matter wave (or resonance), whose existence destroys 
the cloaking phenomenon and makes the 'cloaked' region in fact detectable. However, approxi- 
mate cloaks can be tuned with a precise choice of parameters, close to but not at resonance; the 
flow of the wave from the exterior into the cloak and from the cloaked region out into the exterior 
are balanced, and the cloaking effect is not destroyed, but rather greatly improved [[T3l . We point 
out similar but surface-plasmon based effects [|24l and other subwavelength plasmonics related to 
sensing [I251|26l. 

An approximate QM cloak can be implemented as follows, starting from the ideal 3D spher- 





FIG. 2: The cloak - resonance - Schrodinger hat continuum. Scattering by potentials with different 
parameters, demonstrating three modes: cloak-, resonance- and Schrodinger hat-mode. The graphs show 
the real parts of the effective fields ^^-^-^ on < r < 3. The same curves are shown on different vertical 
scales: (Left) in [-3,3], illustrating the waves outside the cloak, and (right) in [-10,70], where the blow up 
inside the cloak can be seen. (Red) Quantum cloak (with parameter n = rf), for which incident wave 
does not penetrate the cloaked region. (Blue) Almost trapped wave (ri = r[^'^). The cloaking effect is 
destroyed due to the strong resonance inside the cloak. (Black) Schrodinger hat (ri = rf^); probability 
mass is almost entirely captured by the cloaked region, yet scattering is negligible. The red and black curves 
are very close to the incident wave on r > 2, since both a cloak and a SH produce negligible scattering, 
while a resonant cloak is detectable in the far field. The discrepancy shown by the blue curve is due to 
the fact that resonances destroy cloaking and can be observed in the far field; with the parameters here, the 
difference in ^'^•^■^ is small, but the destructive effect on cloaking can be much stronger. For parameters 
used, see discussion of numerical simulations in the Supplement. 



leal transformation optics ElVI invisibility cloak dH. This is based on the 'blowing up a point' 
coordinate transformation [E El y ^-^ x, 



X := F(y) = y, for 2 < |y| < 3; F(y) = 1 



|y| 

2 



forO < |y| <2. 



(1) 



This works equally well in acoustics, forming a cloak with a spherically symmetric singular 
anisotropic mass density and singular bulk modulus [|71 [8l |2I1. Consider the case where the 
anisotropic mass density, M(x), is the identity matrix, and the inverse of the bulk modulus, k(x), 
is = 1 outside the layer 1 < r < 2; the cloaked region is the ball Bi of radius 1 centered at 




ii^ 




FIG. 3: Quantum Three Card Monte. The field ^^™ (left) and the effective field V|// (center) in the 
plane z = when particles are confined in a ball Bl. The central concentration of the wave function 
changes the conditional probabilities of the particles being found in regions in the vicinity of the cloak. 
(Right) The non-normalized probability densities l^'^'^P (red) and iV'ef f P (black) on the positive x-axis. 



origin. For an arbitrary choice of i? > 1, the ideal cloak is then approximated, replacing both the 
mass density and bulk modulus by 1 in the shell (or layer) 1 < r < i?. This gives a non-singular 
mass density Mr and non-singular bulk modulus k^, which approach the ideal cloak parame- 
ters as i? \ 1. Via homogenization theory, the anisotropic mass density Mr is approximable by 
isotropic mass densities m^, consisting of shells of thickness e having alternating large and small 
densities, yielding a family of approximate cloaks [|23l . One then obtains a QM cloak by apply- 



ing the Liouville-gauge transformation ^(x) 



nit 



-1/2, 



x)m(x), so that the Helmholtz equation 



becomes the time independent Schrodinger equation, (—V- V + K — E)il) = 0, where E = cu'^ is 
the energy and Vc = (1 — rrisKjijE + me V- V(m7 ) is the cloaking potential for the energy 
level E. 

For acoustic or EM cloaks constructed using positive index materials, resonances can allow 
large amounts of energy to be stored inside the 'cloaked' region, but at the price of destroying the 
cloaking effect [|22ll23l . However, inserting negative index materials within the cloaked region al- 
lows for the cloaked storage of arbitrarily large amounts of energy; for simplicity, we describe this 
primarily in the context of QM cloaking, where the analogous effect is concentration of probability 
mass. When the cloaking potential is augmented by an internal potential consisting of a series 
of A^ shells, alternating positive barriers and negative wells with appropriately chosen parameters, 
the probability of the particle being inside the cloaked region can be made as close to 1 as desired. 
More precisely, insert into Bi a piecewise constant potential (5(x), consisting of two shells, with 



values Ti, r2, in r < si, si < r < S2 < 1, resp., and zero elsewhere. For suitable parame- 
ters R, e, Sj and Tj of the potential W , we obtain, in the Supplement, a Schrodinger hat potential, 
VsH = Vc + Q. Matter waves incident on the SH are modeled by Schrodinger's equation, 

(-V^ + VsH - E)i, = 0. (2) 

The key feature of Vsh is that the matter waves governed by (|2]) can be made to concentrate 
inside the cloaked region as much as desired, while nevertheless maintaining the cloaking effect, 
quantified as follows. Assume that we have two balls of radius L > 2, one (Bf^) empty space and 
another (Bf'^) containing a Schrodinger hat. Let B'^, B2 denote the central balls of radii 1 and 2, 
resp., for ■ = em or sh, and assume that matter waves ^^"^ and ^^'^ on Bf^, Bf^\ resp., have the 
same boundary values on the sphere of radius L, corresponding to identical incident waves. Define 
the strength of the Schrodinger hat to be the ratio 
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6= ^^ / |^^Vx)|2cix. 

where ^'^'"(x) and ^'*''(x) are solutions which coincide in |x| > 2. We show that, by appropriate 
choice of the design parameters, (3 may be made to take any prescribed positive value. For large 
values of ©, the probability mass of ij)^^' is almost completely concentrated in the cloaked region. 
For £ > 0, the wave ^^^ is rapidly oscillating, and so we also consider the effective wave, V'e//' 
which is obtained as the limit of £ — )• (in a suitable weak sense discussed in the Supplement). 
This is distinct from the 'mirage effect' for standard cloaks, which makes a source within the 
cloaking layer appear to be in a different position due to the chain rule [|20ll27l . 

We next describe some remarkable properties of Schrodinger's hat. To start with, the highly 
concentrated part of the wave function which the SH and the incident wave produce inside the 
cloaked region can be considered as a quasiparticle, which we call a quasmon. A quasmon has 
a well-defined electric charge and variance of momentum depending on the parameters of Vsh- 
Secondly, the amplification and concentration of a matter wave in the cloaked region can be used 
to create probabilistic illusions. Consider (non-normalized) wave functions i})^^ and iIj^^' which 
coincide in Bi — B2, i.e., exterior to the cloaking structure. Then for any region Rout in B^ — B2 
the conditional probability that the particle is observed to be in Rout, given that it is observed 
in Bl — B2, is the same for -ip^"^ and ■ip^'^. However, by choosing the parameters of the SH 
appropriately (see the Supplement), the probability that the particle is in the cloaked region Bi 
can be made as close to 1 as wished. Roughly speaking, the particle tjj^'^'- is like a trapped ghost 



of the particle ip^"^ in that it is located in the exterior of the cloaking structure with far lower 
probability than ^^"^ is, but when ^'^'^ is observed in B^ — B2, all measurements coincide with 
those of V^^"^. When a particle is close to a SH, with a large probability it is grasped by the hat 
and bound into storage within the cloaked region as a quasmon. This is nevertheless consistent 
with the uncertainty principle: although the particle is spatially localized within Bi, the expected 
value of the magnitude of its momentum is large, due to the large gradient of tp'^^' on a spherical 
shell about the central peak; cf. Figs. 1 & 3(center, right). 

The Schrodinger hat produces vanishingly small changes in the matter wave outside of the 
cloak, while simultaneously making the particle concentrate inside the cloaked region. Thus, if the 
matter wave is charged, it may couple via Coulomb interaction with other particles or measurement 
devices external to the cloak. When an incident field ip^" is scattered by the SH, the wave field 
is not perturbed outside of the support of the hat potential Vsh', there are no changes in scattering 
measurements. However, the cloak concentrates the charge inside the cloaked region, proportional 
to the square of the modulus of the value ^*"(0) which the incident field would have had at 
the center O of Bf^ in the absence of the SH. Due to the long range nature of the Coulomb 
potential, this charge causes an electric field which may be strong even far away from the SH. 
If one measures the electric field and the result is zero, then this indicates that ^*"(0) = 0; 
without disturbing the field, one determines whether the incident field vanishes at a given point. 
A measuring device within a Schrodinger hat thus acts as a non-interacting sensor, detecting the 
nodal curves or surfaces on which the incident matter field ^jJ^'^ vanishes, an effect analogous to 
cloaked acoustic and EM sensors ||T3ll24l and near- field scanning optical microscopes [|28l . As 
described in the Supplement, a Schrodinger hat potential also amplifies the interaction between 
two charged particles. 

The behavior of Schrodinger hats and quasmons can be illustrated by means of a quantum 
variant of Three Card Monte, the classic game of chance in which a coin is hidden under one of 
three bowls and the player guesses where the coin is. Consider first a preliminary version of the 
game, played by Alice, who runs the game, and Bob, who makes the guesses. In place of bowls, 
they play the game using A^ empty balls, each a copy of Bf^, and the coin is replaced by a QM 
particle. The surface of each ball is made of material representing an infinite potential wall, so 
that a particle within cannot escape; this corresponds to the Dirichlet boundary condition on the 
boundary. In the game, Alice inserts one particle into one of the balls, after which she mixes the 
balls randomly and asks Bob to guess in which ball the particle is. Bob chooses one and makes 
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internal measurements near the boundary of Bl. Bob, wishing to determine whether the particle 
is in a region R, R C B^ — B2, measures the value of an observable A, which is 1 if the particle 
is observed (X E R) and otherwise; the value of A is also if the particle is not in the chosen 
ball. The expected value of A is p- /i^*" + (1 — p) -0 = p- /i^. Here, p = 1/N is the probability 
that Bob chose the ball that into which Alice inserted the particle, and yU^" = ciem/cem, where 
Oem = /fj |V^^™(x)pc/x, Cgm = || V^'^'"|li2(^ ), and if)^^ is the empty space wave function on Bf^. 

To make the game more interesting, Alice and Bob make a wager: they agree that Bob will pay 
p/x^ € to Alice in advance of each turn, but if he then observes a particle will receive 1 € back 
from Alice. With these rules, the game is fair, with expected profit for both Alice and Bob. 

Now suppose that, before play commences, and unbeknownst to Alice and Bob, a third player 
(the Cloaker) replaces each of the N empty balls with a ball Bf^ equipped with a Schrodinger 
hat. The expectation of A is now p ■ ji^J^ + (1 — p) • = p • fj,^^, where yU^ = ash/csh, where 
ttsh = jp^ \i>'^^{^)\^d-^, Csh = \\'^^'^\\'l2/b )' ^^^ '^^^ ^^ the wave function on Bf". 

Since a Schrodinger hat is an effective cloak, ^^"^ = ip^'^ outside of the ball B2 which contains 
R, and so aem = CLsh- On the other hand, the presence of the Schrodinger hat amplies the wave 
function in Bf'^ and so Csh » Cem', hence ^^^ « /i^^, cf. Fig. 3. When the game is played many 
times, Bob's expected chance of observing the particle in the ball which he chose is smaller than it 
was before the Schrodinger hats were inserted. In other words, after the Cloaker's intervention the 
particles start to disappear from Bob's observations and Bob starts to lose; Alice is unknowingly 
'cheating'. The game can be made as unfair as one wishes by choosing parameters so that 6 is 
very large, using general Schrodinger hat potentials as described in the Supplement. 

We conclude by describing one possible path, discussed in more detail at the end of the pa- 
per, towards a solid state realization of a quantum Schrodinger hat, utilizing a sufficiently large 
heterostructure of semiconducting materials. By homogenization theory, the SH potential can be 
approximated using layered potential well shells of depth —V- and wall shells of height V+. By 
rescaling the x coordinate we can make the values V± smaller (note that in such scaling the size of 
the support of the SH potential grows and E becomes smaller). This sequence of spherical poten- 
tial walls and wells can be implemented using a heterostructure of semiconducting materials. In 
such a structure the wave functions of electrons with energy close to the bottom of the conduction 
bands can be approximated using Bastard's envelope function method [|29ll . Choosing the materi- 
als and thickness of the spherical layers suitably, the envelope functions then satisfy a Schrodinger 
equation whose solutions are close to those corresponding to the SH potential. 

8 



Supplemental Material 



In this supplement, we provide the rigorous analysis needed to confirm the existence and 
behavior of Schrodinger hats, specify the parameters used in the Figures, and detail the 
proposed solid-state implementation. 

ANALYSIS OF QUANTUM CLOAKS AND SCHRODINGER HAT POTENTIALS 

An approximate acoustic cloak 

Below we will use the approximate cloaks modeled by the Helmholtz (or the Schrodinger) type 
equation 

(-V- M^ V - u^kr) ur = in the domain Q C ]R^ (3) 

ur = h on the boundary dil, 

where J? > 1 is a parameter corresponding to the effectiveness of the cloak, to is the frequency, 
and Mr and kr are the coefficient functions defined below. Let u denote the outward unit normal 
vector of dD.. Measurements on the boundary dil are mathematically modeled by the Dirichelet- 
to-Neumann operator defined by 



AR{h) = u ■ M^'VurI 



dn, 



describing the response of the system, i.e., the Neumann boundary value u ■ M^^VurIqu, when 
the Dirichlet data ur = h is posed on the boundary. In the theory of the approximate cloaks the 
coefficient functions Mr and kr are constructed in such a way that as i? \ 1 the Dirichlet-to- 
Neumann operators A^ approach to the Dirichlet- to-Neumann operator A^°™°^ : /i i— )■ z/ ■ Vwlan 
for the boundary value problem 

(-V- V - w^) M = in the domain fi C ]R^ (4) 

ur = h on the boundary dVt, 

modeling empty space. In practical terms, this means that when the parameter R is close to 1, 
for the approximative cloak all boundary observations on dVt are close to the observations on dVL 
made when the domain Vt is filled with a homogeneous, isotropic medium. 



Approximate cloaks are the basis of our construction of Schrodinger hat potentials. We start by 
recalling some facts concerning nonsingular approximations to ideal 3D spherical cloaks [|22l l23l 
[30] - l33l . For R > 0, let Bji = {|x| < R} and Sr = {|x| = R} be the open ball and sphere, resp., 
centered at the origin O and of radius R in three-space. Moreover, let Br = {|x| < R} be the 
closed ball. For 1 < i? < 2, set p = 2(R — 1), < p < 2, so that i? \ 1 as p \ 0, and introduce 
the coordinate transformation Fr : Bl — Bp —^ Bl — Br, 

{y, for 2 < lyl < L, 

(l + Mj^, forp<|y|<2. 

For i? = 1 (p = 0), this is the singular transformation of [[Tl|2l|4l, leading to the ideal transforma- 
tion optics cloak, while for R> 1 (p > 0), Fr is nonsingular and leads to a class of approximate 
cloaks [|22l l23l [30l [32l [33l . Thus, if Mq = 5jk denotes the homogeneous, isotropic mass ten- 
sor tensor, then, for i? = 1, the transformed tensor becomes an anisotropic singular mass tensor, 
Mijfc(x), on 1 < |x| < L, defined in terms of its inverse. 
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-l^i'^.„^ ..^^ .^,;^.,.^ ^ >A 9(Fi)\„, 9(Fi 



\k 



■ (6) 

x=Fr'(y) 



This means that in the Cartesian coordinates Mi (x) is the matrix with elements 

(Mi),fc(x) = ^{6,, - P,,(x)) + ^(|x| - l)-2p,fc(x), K |x| < 2, 

where the matrix P(x), having elements Pjk{x) = \x.\~'^XjXk, is the projection to the radial direc- 
tion. 

On the other hand, when R > 1, we obtain an anisotropic but nonsingular mass tensor, 
MHjfc(x), on i? < |x| < L, given by 



(M^y (x) = {{FR)jy'iy) := i V ^^(x) ^^(x)<5^^(x) 

\ R J K J yy nj* J u; det[^(x)l ^ dx^ ^ ' dx'i ^ ' ^ ' 

^ ax ^ 'J p,q=l 



•(7) 



For each R> 1, the eigenvalues of Mr are bounded from above and below; however, two of them 
/^ CX3 as _R \ 1. We define an approximate mass tensor Mr everywhere on Bl by extending it as 
an identity matrix, 

, Mr.M fori?< Ixl < 2, 

5jk, for |x| < i? or 2 < |x| < L. 
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In sequel, we use the notation M^j .,fc(x) also for M^^L(x). We define a scalar function /€i?(x) on 



i^rW 



ri{x.), for |x| < i?o, 

64|x|-^(|x| - 1)^ for i? < |x| < 2, (9) 

1, for i?o < |x| < i? or 2 < Ixl < L. 



where 

N 

77(x) = ?7(x;r) = ^TjX{s,_r,s,){\^\). (10) 

i=i 

Here r = (ri, T2, ■ ■ ■ , t^), Tj G M are parameters which one can vary, = sq < Sj < s^ = Rq 

are some fixed numbers, and X(s _i,s )(^) is the indicator function of the interval (sj-i, Sj). This 

means that we have a homogeneous ball Bg^^ coated with homogeneous shells. Sometimes we 

denote /€r(x) = /€i?(x; r). Note that in acoustics kr has the meaning of inverse of bulk modulus; 

later, in quantum mechanics, it gives rise to the potential. 

Below, we consider what happens as i? \ 1. In fact, for rigorous mathematical analysis we 
should modify the above definition of kji by replacing it, e.g., by ki so that for all x E B^ the 
quadratic form corresponding to operator J^ (Vf ■ M^^Vv — uj'^Kji\v\'^)d'x becomes smaller (for 
any fixed v) as R decreases. However, in order to compute solutions explicitly and to present 
considerations in a simplified way, we will consider the case when kr is defined as above. Math- 
ematical proofs will be presented elsewhere. 

Next, consider in the domain Bl the solutions of the Dirichlet problem, 

(-V- M^ V - w^^r) Mr = in Bl, ur\s^ = h. (11) 

Since, for 1 < i? < 2, the matrix Mr is nonsingular everywhere, across the internal interface Sr 
we have the standard transmission conditions. 



ur\sj,+ =Ur\s^_, (12) 

er- {M^'Vun) \s,^ = e,- [M^'Vur] \s^_, 

where er is the radial unit vector and ± indicates the trace on S^as r ^ R"^. 
In the physical space B^ one has 



;+(F^i(x)), fori?<|x|<L, 



uni^)={ -^ - — ' ' (13) 

f^(x), for |x| < R, 
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with v^ in the virtual space, which consists of the disjoint union {Bl — Bp) U Br, satisfying 

(_V2 - u^)v+{y) = for p < |y| < L, 



and 



-V^-tu\R)v^iy)=0, for|y|<i?. 



(14) 



With respect to spherical coordinates (r, 6, (p), the transmission conditions ( 12 1 become 



p^drV^{p,e,<f))=R^drV],{R,e, 



(15) 



Since Mr, kr are spherically symmetric, cf. ( |8|9p , we can separate variables in ( [TT] ), representing 

ur as 



oo n 



UR{r,e,<P) = J2 E <'"W^r 



'^'"^-^-7(^,0), 



(16) 



n=0 m,=— n 



where Y^ are the standard spherical harmonics. Then equations (11) give rise to a family of 



boundary value problems for the u^™. For our purposes, the most important one is the lowest 
harmonic term (the s-mode), Uj^ , i.e., the radial component of ur, which is independent of {6, 0). 
This is studied in the next section. 

Spherical harmonic coefficients 



The lowest harmonics. For i?o < 1 < -R < 2, consider the Dirichlet problem on the ball Bl, 



{-V-M-'V-uj^kr)ur = in Bl, 
ur\sl = ^(x). 

We will express asymptotics in terms of the quantity p = 2(i? — 1) as p \ 0. 
We have shown elsewhere [|23l that 



(17) 



For a specific value of the parameter r E M^, denoted r = t'^'^^{R, co), there is a blow-up 
effect, or interior resonance, destroying cloaking. This corresponds to the case when r is 



such that there equation ( frTj ) has a non-zero radial solution with h = 0. In this case the 
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solution ur grows very much inside the cloaked region as i? — )■ 1. his means that the inside 
of the cloak is in resonance and the wave tunnels outwards through the cloak, so that this 
resonance is detected by boundary measurements outside of the cloak. 

For another specific value of r, denoted r = r^'^"(i?, to), the cloak acts as an approximate 
cloak and inside the cloaked region the solution is proportional to the value which the field 
in the empty space would have at the origin. This corresponds to the case when the equation 



( 17) has a radial solution ur which satisfies drUpt{L) = ujj'^iujL) and uj^{L) = Jq{luL), or 



equivalently, ur{x) = jo{co\x.\) for R < |x| < L. 



Due to the transmission condition ( 12) we see that the values t'^^^{R, u) and t^^'^{R, uj) are close 
and lim/j-j.! t^^"-{R, u) = limij_5.i r^'^^{R, u). 

We now explain in detail how to choose r = t^'^'^{R,uj): First, fix R and to, and choose 
< i?o < 1- Consider the ordinary differential equation corresponding to the radial solutions u(r) 



of the equation ( 17 ), that is. 



- ;:^^r'^«W^"(^)j -^^f^RirMr) = 0, (18) 

and pose the Cauchy data (i.e. initial data) at r = L, u{L) = Jq{Luo), dru{L) = uJq{Lu). Here, 
crij(r) is the rr-component of the matrix Mj^^, that is, cr/j(r) = 2(r — 1)^, for i? < r < 2 and 
<^Rii^) = !> elsewhere. Then we solve the initial value problem for the the ordinary differential 



equation ( 18) on interval r G [Rq, L] and find the Cauchy data {u{Ro),^{Ro)) at r = Rq. Note 



that on the interval r G [Rq, L], i^r does not depend on r. Consider next the case when 

AT = 2, si=Ro/2, S2 = Ro, r = (ri,r2), (19) 

where ri and T2 are constructed in the following way: 

First, we choose T2 to be a negative number with a large absolute value. We then solve of the 



initial value problem for ( 18) on interval r G [si, Ro] with initial data {u{Ro), ^{Rq)) atr = Rq 



In particular, this determines the Cauchy data (m(si), ^(si)) at r = Si. 

Secondly, consider ri, r2, as well as R, Rq, to be parameters, and solve the initial value problem 



for (18) on interval r G [0, si] with initial data (m(si), 4^(si)) at r = Si. Denote the solution by 



u{r; Ti, R, Ro, T2) and find the value |^(r; ri, R, Rq, T2)\r=o- Then, for given R, Rq, and r2, we 

find Ti > satisfying 

nil 

— {0-n,R,Ro,T2)=0. (20) 

ar 
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We choose ti to be the smallest value for which (20) holds, and denote this solution by 
Ti(R, Rq, T2). Summarizing the above computations, we have obtained a cloak at the frequency 
CO, that is, for the energy Eq = co"^, such that its radial solution m(x) satisfies m(x) = jo(i^|x|) 
for |x| G [2, L]. Moreover, when r2 is large, this solution M/j(r) grows exponentially fast on the 
interval [si, S2], as r becomes smaller, while on the interval [0, si] it satisfies ^^(0) = 0, so that 



u{r) defines a smooth spherically symmetric solution of ( 17 1. 

In the context of QM cloaks below, the construction above can be considered as follows: Inside 
the cloak there is a potential well of the depth ri, enclosed by a potential wall having the height 
r2. The parameters ti and r2 are chosen so that the solution is large inside the cloak due to the 
resonance there. Moreover, the choice of the parameters is such that the flow associated to the 
wave function from the outside into the cloaked region and from the cloaked region to the outside 
are in balance. The cloaked region is thus well-hidden even though the solution may be very large 
inside the cloaked region. In a scattering experiment, with high probability the potential captures 
the incoming particle, but due to the chosen parameters of the cloak, external measurements cannot 
detect this. 

Using the implicit function theorem, one can show that for generic values of T2 and 
Rq, there is a limit UmR\^i ri(i?, i?o, T2) = ti{Rq,T2). We note that the solution uji(r) = 



u{r;Ti{R,Ro,T2),R,T2) of (IS) has \imit\imji\iu{r;Ti{R,Ro,T2),R,T2)) = c^{r), r < 1, 



where c G C and $(r) ^0 an eigenfunction of the boundary value problem 

(V^ + cu^K.ir, r))<l>(y) = 0, for |y | < 1, 

a,$(y)|,=i = 0, (21) 

where r = (ri(l, Rq, T2), T2) and $ is normalized so that ||$||l2(^^) = 1. 

Higher order harmonics. Let r = {ti(Ro,T2),T2). As ti was chosen to be the smallest 



solution of ( 20 1 we have that uj an eigenfrequency of the problem 

(_V2-a;V(r;r))t;(y)=0, for |y| < 1, 



9,t;(y)|,=i = 0. (22) 



Let us now analyze the solution ( 13 ) using spherical harmonics. Recall that in B^ the function 



f^(y) = Mij(y) is a solution to the homogeneous equation ( 14). Thus, in particular, 



'^Rir,9,^) = J211 uT{r)Y:\9,y,), (23) 



n=0 111=— n 
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where 



R (r) = anmjn{u;r) +pnmh^n\^r), Rq < r < R, (24) 



wTir 



M^"^(r) = a^rnJniujVrir), < r < Ro/2, 

and ^/T2 is pure imaginary. Here a^^ = anm(w; -R), ctnm = anm(w; R), Pnm = Pnm{^] R) etc. 
are yet undefined coefficients. Note that the terms with hniuy/r^ r) are absent near r = since 
t>^(y) has no singularity at 0. 
Now, for p < r < L, 



00 n 



v^[r, 



^'^) = E E {Cnmh';^\ujr) +bnn.Jn{cOr))Y:^{e,y^), 



n=0 m=— n 

with as yet unspecified bnm = hnm{'^] R) and c„m = Cnm{'^] R)- 

Expand the boundary value h on dBj^ in surface spherical harmonics as 



00 n 



h{d,^)=J2 E fnmYr{0,V). (25) 

n=0 rn=—n 

As shown in the previous section, 600 = foo, cqo = 0, and both poo and aoo can be solved for using 
the transmission conditions, which determines the coefficients for n = 0. 

Next we consider the higher-order coefficients, for n > 1. To simplify notations, denote by 



U^"^{r)Y;^{e, ifi) the solution of (|T4|) for which [/^""(r) = jn{(^^/nr) for r < Ro/2. Then the 
coefficients in (23) can be written in the form u^{r) = anmU^{r). Observe that there exist the 



limit 

lim U'pi.r) = [/{""(r) 
and due to the way the coefficient ti = ti{Rq, T2) was chosen, for generic values of T2 we have 

drU'^'^il) ^ for all n > 1, -n < m < n. (26) 



Next we assume that T2 is such that (26 1 holds. 



By the transmission condition ( 12 1 



ar,m{R)Ur{R) = bnmiR)JnicOp) + Cn^{R)h^^\uj p) , (11) 

R'a^m{R)drUr{r)\r=R = pXm{R)drJn{uJr)\r=p + Cnm{R)drh^n\p)\r- 
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=p- 



Using asymptotics of Bessel and Hankel functions [|34l . we obtain from (27 1 that 

Cnm{R) = — — — — ya„„(i?)5,f/r(l)^"+V" + 0(p"+^)), (28) 

[2n + l)[2n)\ 



Note that above drUi^{l) is non-vanishing by (26 1. Using (28 1 and the Dirichlet condition (25 1, 
we finally see that 

a^m{R) = (2n + l)«„/„„a;'^— ^p"+i + 0{p^^^). (29) 

Together with the transmission conditions onr = Rq and r = Rq/2, this implies that 

h^n = -^r^fmn + 0(p), Cmn = 0(p'"+^), (30) 

(^mni (^nmi (^nm ^\P )) Pmri) Pnm ^\P )• 

The above considerations for n = and forn > 1 can be summarized as follows: As i? \ 1, 
the solutions m^(x) converge in Bl — -B2 to the solution u corresponding to the homogeneous 
virtual space, 

{-V-V-uj^)u = in Bl (31) 

and, in the domain Bi, to the solution in empty space, 

lim u^(x) = /3u(0)$(x) (32) 



where $(y) = $(r) is the radial solution of the equation (22), /3 = -^ and u{0) is the value of 



*(i) 



the solution of (31 1 at the origin. 



Next we consider the implications of this for quantum mechanics. 

Approximate isotropic cloaks in quantum mechanics - Schrodinger's hat potential 

The approximate anisotropic cloak (Mr, kji) can be further approximated by an isotropic cloak 
{rriR^e, f^B), where m/je(x) is a smooth isotropic (i.e., scalar-valued) mass density, which we de- 
note by lowercase m to distinguish it from the anisotropic mass tensor denoted by M. It satisfies 
Ci(i?, e) > rriji^si'^) > Ci > 0, and leads to an approximate cloak equation, 

^ -V - u^kr] ur. = mBL (33) 



mR^e / 
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We will use isotropic mass densities which, for J? < |x| < 2, are of the form 

1 ,, ,, /Ixl— i?\ ,,, ,, /Ixl— i?\ f\x.\—R\ ^^^^ 

a{\x\)pi ^ +&(|x|)p2 — +P3 — . (34) 



mR^sK^} \ e J \ e J \ ^ 

Here, pi, p2, Pa are bounded non- negative smooth functions with period one such that pi, P2 = 
near integer values, while ps = 1 near integer values. For each p > we choose a sequence of 

^k = £k{p) — )• as A; — 7- oo such that |y^ is an integer. As for |x| < R and |x| > 2, we take 
mp£(x) = 1. It is possible to choose a(|x|) = 0(1) and fe(|x|) = 0(|x| — 1) as |a;| \ 1, so that 

1. The mp£^.(x) are smooth functions in Bl\ 

2. The nip^eki'^) approximate Mi as er^ — )■ and then p — )■ 0. Namely, the operators 
— V-m^^V converge (as described below) to — V- Af^f^V as e^ — )• and p — )• 0. Note 
that these a(|x|), &(|x|) are independent of p. 

Below we use the shorthand notation e — )■ instead of Sk -^ 0. Denote 

^(x) =a(|x|) f pi(r')rfr' + 6(|x|) f p2{r')dr' + f p^{r')dr\ K |x| < 2, (35) 

JQ Jo Jo 

^(x) = 1, |x| < 1 or |x| > 2. 

We can choose Pi{r'), i = 1, 2 to be smooth functions that are very close to the characteristic 
functions of the intervals (0, 1/2) and (1/2, 1), continued periodically, while p3(r') has its support 
very close to r' = 0, continued periodically. To have that the Hamiltonians corresponding to irip^^ 
to approximate Hamiltonian corresponding to Mi (cf. [|23l for analysis of the F-convergence and 
the two- scale convergence of these operators) we need that 

(a(r) + 6(r)) ^ 2, - ^^ + 



2' ' ' ' '' '2 \a{r) b{r) J 2(r-l)' 
i.e., 

a{r) ^ 2(1 + V2-r), b{r) ^ 2(1 - V2-r), 1 < r < 2. (36) 

Thus, the mass density mp e(|x|) corresponds to two materials occupying layers of equal width in 
each cell with those cells separated by a very thin layer of uniform density. 



Making a Liouville gauge transformation for equation ( 17 ), that is, introducing 

V'p,e(x) =m;y^(x)u/j,£(x), (37) 
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(33) becomes the Schrodinger equation with potential V^^, supp(V^e) C B2 — Bi. Here, 
supp(Vp,£), the support of Vp^e, is the set where Vp,e(x) is non-zero. For generic values of the 
parameters, these model approximate invisibility cloaks for matter waves [l22ll . with Schrodinger 



equations 

(-V ■ V + ^p,, + Qp - E)Vp,. = mBL, u\dB^=h. (38) 

Here E = w^, and these cloaking potentials, V^^^, are of the form, 

V^p,,(x) -E = mJ/jV ■ V(m^'f) - EmR^.nn, i? < |x| < 2, (39) 

Vp^ei^) =0, < |x| < i? and 2 < |x| < 3. 



In addition, inside the cloak, i.e. in Bi, there is the cloaked potential (cf. ( 10)) 



gp(x) = -^(r/(x; r) - 1), supp(Qp) C 5^^ C B,, (40) 

where r = t^^^{R, co) and Qp depends on the parameter t(R, cu). Thus, when r is appropriately 
chosen, the total potential Vp^^+Qp in B^ acts as a quantum mechanical sensor cloak. It is this total 
potential which we call a Schrodinger hat (SH). (The cloaked potential, that is, Qp, resembles a 
"hat", specifically the Mexican hat, but the terminology is chosen because of the interesting effects 
of the potential, not because of its profile.) 

Analysis of Schrodinger's hat: convergence of normalization constants 

To consider the properties of a Schrodinger's hat, one first needs some convergence properties 
of approximate cloaks as they approach an ideal one, initially treating the case when E is not a 
Dirichlet eigenvalue of Laplacian in the ball of radius L. The case of Dirichlet eigenvalues will be 
considered later. 

Let u satisfy the boundary value problem 

{-V-V -E)u = 0, in5LCM^ (41) 

uldBr = h. 



Moreover, let 



Qq = limQp = E(l -r](^;(Ti(Ro,T2),T2))), 
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(see (40)), be the potentials supported in S/?,,, with the parameters suppressed. Then the Qo 
describe QM cloaked sensor potentials at energy E. 

Let F = Fji with R = Ihe the singular blow up map and define 

m(F-1(x)), forxG^L-^i, 
/3m(0)<I>(x), forxG fii, 

where /3 = $(1) and u is the solution of (41 ). 



MX 



(42) 



Let Up be a solution of Helmholtz equation 

UpldBL = h, 

and Up £ be a solution of Helmholtz equation, 



0, in Bl C 



(43) 



(- V- m;jV + Qp- EKR)up,e = 0, in Bl C 
Up,£\aBL = h. 



(44) 



Then, 



limti 

e-s>0 



p,e 



U 



pi 



(45) 



in the space L'^{Bi). Moreover, letting K C B^ — Bihc arbitrary, using (35 1, (37 1, one can show 
that 



limlini / |'?/'pe(x)P(ix = limlim / ?7ipe(x) '^\ups('x)\'^dx= / 9hi.)\uCx)\'^dx. (46) 



Proofs of the claims (45) and (46) are omitted and will be detailed elsewhere, but the principal 
ideas are that, using the above properties of the spherical harmonics expansions. 



limlimM„e(x) = m(x) for x G B^ — dBi, 

p-s>0e-5-0 



and that, for G C{Bl), 



p— >0£— 5>0 J j^ 



limlim / r7i^(x)0(x)(ix 



0(x)0(x) dx. 



(47) 



(48) 



K 



In the following, we use 



^' 



effi 



0(X 



,1/2 



MX 



(49) 



which can be considered as the effective wave function, when one is modeling the location of the 
particle corresponding to the wave function Wp ^ with sufficiently small p and e. 
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Conditional probabilities for location of particles in a Schrodinger hat 

Consider next the situation of two balls of radius L, one empty and the other with a Schrodinger 
hat potential, and assume that there is one particle in each ball. Let n(x) be wave function corre- 
sponding to the particle in the empty ball and ^p e be the wave function corresponding to the ball 
with the SH potential. Assume that u and ipp^s have the same boundary value, h, on dB^. 

Consider next the normalization constants 

Jbl-B2 

cr^'y= [ _ iiz(x)|2rfx, 

Jb2-Bi 

Cf""^ = lim lim / iV^p .(x) |2 dx, 

Cf '^'^ = lim lim / |^p,.(x)|2cix, 

Cf -"^^ = lim lim / |^p,.(x)|2rfx. 

The probability that a free particle X in the empty ball B^, described by the wave function u, is in 
fact located in the smaller ball Bi, is equal to 

^ empty 
\\ 1// y^empty . y^empty . /^ empty ' 

Similarly, the probability that a particle Xp ^ in the ball Bl with the SH potential V^ ^ + Qp, 
described by the wave function i/jp^^, to be in Bi satisfies 

/^cloak 

limlimP(|X„, G BA) = — p^ V^ tt^^- 

Similar results hold for the probabilities for the particles to be in B2 — Bi and Bl\ B2. 



Using (46)-(48), we see that 



ricloak ^empty 

Recall that u is the solution of (41 1 and $ is the L^(i?i) -normalized Neumann eigenfunction inside 



the cloak, cf. (42). As for Q'°"'=, by (46) and (49), we have 



Cf'^'= [ \r"{^)\'d^= f ^(x)|M(x)|2t;x= / 9{y)\u{y)\'dy, (50) 

JBoXB, JBnXB, JB2 



'B2\Bi JB2\Bi JB2 
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^(y) = ^(x) ^ , x = F(y), |y| < 2. 



From the definition of F{y), cf. (5 1, one sees that ^(y) = O ( pW ) ; however, due to the bound- 



edness of the norm of u in the space H^{Bi), the integral in (50) is bounded. An important 
observation is that, by changing the parameters r determining the potentials Qp, Qo, one can keep 
Cf°"-^ and Cf""'^ unchanged while Cf""-'' is made arbitrarily large. 

Consider a particle X with energy E in an empty ball Bl and another particle Xp^ with energy 
E in the same box but with the SH potential Vpe + Qp. Let us next consider the event that the 
particle X in the empty ball corresponding to wave function u is located in the set S <Z Bl, and 
denote this event by Ls- Then, 

^\^S) = Cempty / p(Xj| GX, Cempty = „empty , ^empty . ^empty ' 

The conditional probability for the event that the particle X is in 5 C Bl\B2, conditioned on the 
particle being in So := Bl\B2, is thus 

(' \u(x)\'^ d:si 

nLsLs,) = ;^ ) '', ■ (51) 

Next, for a ball B^ in which there is an SH potential, denote the event that a particle Xp ^ is 
located in the set 5 C i?L \ -B2 by Ls{p, e). Then, 



Js 
limlim Ccioak{p,£) 



p^Oe^O ^'°'^'=^^' ^ (jcloak ^(jcloak ^(Jdoak- 

Thus, the conditional probability for the event that the particle X is in 5 C B^ \ B2, conditioned 
on it being in 6*0 := -Bl \ -B2, is 

p(Z,.(p..)iZ,.>..))^ |i,^;;Wg- . (52) 

and 

hm\imF{Ls{p,e)\Ls,{p,e))=FiLs\Ls,). (53) 

In summary, the above computations have the following consequences for the wave functions u 
and ijjp^e, corresponding to the the particles in an empty ball and a ball with the SH potential, resp. 
Recall that the both field u and ^p^^ have the boundary value h. Then, under the condition that we 
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observe a particle in Bl \ B2, the conditional probability that the particle is observed to be in a set 
S C Bl \ B2, conditioned on it being observed to be in a set Bl\ B2, is same for both balls. 
Consider now the case when S = Bi, or, more generally, S C -Brq C Bi. Clearly, equations 



([5T])-(p2|) remain valid for such S. However, equation ( [53] ) is no longer valid. Moreover, by 
choosing properly parameters r, we can make limp^o lini£-i.o ^{Ls{p, £)\Lso (p, £)) as large as we 
want. To measure this effect we introduce the ratio 



p^os^o\u{0)\^ Jbi |$(1)|^ 



see (42), (49), where 9(x) = 1 in i?i; we call & the strength of the SH potential Qq. & depends 
only on the choice of the parameters t{p, u) = (n, r2), in particular, the parameter T2 that deter- 
mines how rapidly the solution grows in the layer Bs^ \ Bg^ C -Br,,. Choosing t(p, lo) appropri- 
ately, we can achieve any prescribed positive value of &. The parameters (ri, T2) of the SH poten- 
tial do not change the (non-normalized) wave function outside the ball Bi but change radically the 
wave in it. The wave inside the ball Bi has the variance of momentum © |n(0) p J^ | V$(x) prfx 



and the charge Q' considered later in formula (69). We consider the function inside Bi as a quasi 



particle and introduce a solid state model for in Section [} 

We analyzed above the eigenfunction in a ball, and it was enough to analyze only the lowest 
harmonic. One can replace the ball B^ with an arbitrary domain fi C M^ containing the ball B^ 
where the cloaked SH potential is supported using the following observation: The boundary value 
problem 

(-V- V + Vp^e + Qp- E)ipp^^ = = in the domain Q, (55) 

^p,e = h on the boundary dVl, 

is equivalent to the problem 

(-V- V - E)ijp^e = in the domain Q - Bl, (56) 

■ipp,e = h on the boundary dVl, 

'^■'^'^P,e\dBL = KA'^P,e\dBL) on the boundary OBl , 

where Ap ^ is the Dirichlet-to-Neumann operator for the equation (—V- V + Vp^ + Qp — E)-ip = 
in the ball B^. The previous analysis for the higher order harmonics shows that the Dirichlet- 
to-Neumann operators Ap^^ for the ball with the SH potential tend, as p and e tend to zero, to 
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the Dirichlet-to-Neumann operator corresponding to an empty ball. Thus, the above analysis of 
the behavior of the solutions in the ball B^ can be readily generalized to an arbitrary domain 
n. Potentials which, for some incident wave, produce a scattered wave which is zero outside a 
bounded set, are said to have a transmission eigenvalue fS5\. We emphasize that the scattered 
wave caused by the SH potential is approximately zero for all incident fields. 

An illustrative example of this phenomenon via quantum three-card monte is described in the 
body of the paper. 

Scattering from the hat and particle storage 

Now consider the effect of Schrodingers hats on scattering experiments in M'^. In the case of 
free space, the wave function ^ satisfies 

(-V-V-E)* = 0, inM^ (57) 

and we can choose 

^(x) = ^^"(x) = e^'^^^''^), |e| = 1, oj'^ = E, 

so that ^*"(x) are the plane waves in the direction e. Compare them with the wave functions in 
M^ corresponding to scattering from the SH potential, 

(_VV + l^,,, + Qp-^)^,,, = 0, inM^, (58) 

where \E'^^g(x) satisfies Sommerfeld's radiation condition. Note that the SH potential V^^ + Qp 
vanishes outside -82- 

As scattering data for these problems is equivalent to Dirichlet-to-Neumann operators on dBj^ 
1361, we can use previous results to consider scattering from the SH potential. In particular, we 



see using (37), (42), and (45), that when p and e are small enough, the solution \l'pe(x) is close 
to ^^"(x) outside B2, close to m;,£/^(x)^*"(F-i(x)) in B2 - Br, and close to ^*"(0)$(x) in the 
cloaked region. 

Thus, we see that scattering observations, i.e., observables depending on the far field patterns 
of the solutions, are almost the same for the empty space and for the Schrodinger hat potential 
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(when p and e are small enough). However, if we consider the conditional probability of particles 
in an empty ball and one with with the SH potential, i.e., the ratios 

4iv^(x)prfx 



'51,52 



and 

/5jV^p,.(x)prfx 



-^Si,52 



/5J^p,£(x)Po?x' 



we see that one can make Isi,S2 >> ^Si,S2 if ^i contains Bi and 5*2 does not intersect B2. The 
physical interpretation of this is that, when a particle scatters from a Schrodinger hat, there is a 
high concentration of probability mass in the cloaked region, but this can not be detected from far 
field observations. 

A hat and two ions - an interaction amplifier 

A SH potential also amplifies the interaction between two charged particles. To see this, con- 
sider two particles of the same species, which for simplicity we simply refer to as ions. Suppose 
that one ion, having energy E, is confined to a domain which also includes a SH potential Vsh, 
centered at point O. Tuning the parameters of the SH potential, the probability that the ion is 
concentrated in a ball B^, of small radius d and center O, can be made as close to 1 as desired. If 
one now inserts two ions into the domain which interact both with the potential Vsh and with each 
other via a Coulomb potential, both charged particles can not be concentrated in Bd- First-order 
perturbation theory shows [? ] that the energy level of the particles in the presence of a SH behaves 
as if the charges of the ions were multiplied by a (large) factor 0{d^^/'^). 

Let us consider a large ball B^ (we could also consider a box D^ = [— L, L]^ or any other 
domain) containing a SH potential and a single particle system modeled by 

(-V- V + \/p,,(x) + Qp(x) - ^)^p,,(x) = in Bl, iPpMl = 0, UpA^iB,) = 1, (59) 



cf. (38 1. Here .E > is the Dirichlet eigenvalue of the partial differential operator in (41 1 in the 
empty ball B^, corresponding to radially symmetric eigenfunction ^°(x) = Jq{uj\x\), u = \/E, 
so that jo{ujL) = 0, and V^,e(x), Qp{:>i.) are chosen as in the subsection. I.C. In addition to the 



24 



one-particle equation, we consider a system of two charged particles, modeled by 

(-v.- V. - Vy Vy + \/p,,(x) + Qp(x) + y,,,(y) + Qp{y) + ^^—^ - 2E')^{^, y 



x-y| 



in(x,y) eBlX Bl, 



^0, 
(60) 



where the energy 2E' is close to 2E. Here, a is proportional to the product of charges of the ions. 
When a is small (which we now assume), the solution \1' is a perturbation of \1'°, the product of 
two one-particle solutions. 



v[/0(x,y)=^,,,(x)V^,,,(y). 



(61) 



By first-order perturbation theory, we can write 



E' = E + aEi + 0{a^ 



(62) 



v[/(x,y) = vl/W(x,y) + avl/(;)(x,y) + 0(a2), ^«|a(B,xB,) = 0. 

We will show that the two-particle interaction is strongly effected by the presence of a SH. 

Substituting approximation (62) in equation (60) and considering terms of order 0{a), we 
obtain the equation 

( - V,.- V, - V,- V, + i;,,(x) + g,(x) + V,M + Q,(y) - 2E)v|/W(x, y) (63) 
= (^2Ei - ^^) ^i?(x,y), for (x,y) e Bl x Bl, ^i;]|o(B,xs,) = 



,(0) 



As 2E is the eigenvalue of the left hand side of (63 ) with \l'p,i being its eigenfunction, it is neces- 
sary that 



El 



B^ Jb^ |x - y 



#P,£(x)n^p,e(y)P rfxc^y, 



(64) 



where we have used (61 ), (59 ). Let 



$i(x)=/ vl;^,,(x,y)^,,,,(y)rfy, $i'Hx) = / vl/(i)(x,y)^,,,,(y)c?y- 
JBl JBl 



Multiplying (63 ) by ^/^^ ^(y), integrating over i?L in y, and using integration by parts yields 



(-V,- V, + l^p,,(x) + Q,(x) - E) $i^Hx) 
2-Ei^p,e(x) - V;//(x)^p,£(x), 



(65) 
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with the effective potential given by 

K//(x)=/ ^^|^,,,(y)|%. (66) 

Next, we compare the above results with the case when a single particle scatters from the potential 
which is the sum of to the SH potential V^p,e(x) + (5p(x) and the potential Ve//(x) multiplied by 
parameter a, that is the equation 

(-V- V + V,,e{^) + gp(x) + aVeff{^) - E'ff) 0^^^(x) = in Bl, (p'^^lan, = 0. (67) 

First-order perturbation theory then implies that 

E^// = E + 2aEi + 0{a^) = E' + ©(a^), 

r"{^) = ^p,,(x) + a<l>^l\^) + Oia') = <^l{^) + 0{a'). 

Thus, when two particles are in a ball containing a SH potential, each particle behaves, up to error 
0(0^), as if the other particle and the SH potential were replaced by the potential Vp,e(x) +Qp{x) + 

aVeff. 

To analyze Vejf, let & be the strength of the SH potential, as defined in (54). Note that by 
choosing r appropriately, we can make (5 arbitrarily large and, for such (5, the L^ normalized 
waves are strongly concentrated in fi^. For large & and every x E B^, 

f -^\^,M\'dy«f ^^|^p,.(y)|%. (68) 

J BlXBh |x — y| Jbh |x — y| 

Recall that the potential Qp is constructed so that the solution inside the cloak is concentrated 
in a ball Bji^, and assume next that -Rq > is small. To emphasize this, we denote Rq = 5 « 1. 
One then obtains 

K//(x)^ / ^-|^,,,(y')|W^rT' Q'= [ l^pAyTdy' (69) 

Jb, |x-y'| |x| Jb, 

Thus, $l(x) satisfies a one-particle Schrodinger equation where the potential is a SH potential 
slightly modified by a Coulomb one with charge aQ' at the origin. 

Compare this with the case when we have no SH potential but add the Coulomb interaction. 
Analyze this using first-order perturbation theory, writing the wave function and the energy as 

$f^°"(x) = (^<^'"P*^(x) + a(j)^°^(x) + 0(a'^), 
E^'^ = E + aE^'^ + 0{a^), 
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where 0'''"p*^(x) = CLJoiuJo\x\) and cl > is such that ||0''"'^*^||l2(b^) = 1. Analogously to ( [65| ), 
one obtains 



■V ■ V0r"(x) - S^°"0i "(x) + aV(x)0^°"(x) = O(a'), in Bl \ 0, (70) 

Cou\ f\ 

1 |9Bi - U, 



where 



JBr^ |x - y I 



and 



^ JBi JBi |x - y| 

Now compare Ei defined for the Coulomb and SH potentials with the analogous quantity i?f °" 
defined for the Coulomb potential (with no SH potential). Consider the case when & and cl 
are of the same size, and when the potential Q^ is constructed so that the solution inside the 



cloak is concentrated in a ball Bs, where 5 > is small. Then formula (64) implies that we 
have El = 0(6^^). When 6 is very small, the value aEi, i.e., the change in the energy level, is 
much larger than in the Coulomb case without SH potential. Thus, by engineering the potential 
Qp appropriately, one can make the SH cloak increase the interaction caused by the Coulomb 
potential: When two particles are put in a ball with a SH potential, the energy level of the particles 
is changed as if the charges of the particles were multiplied by a factor of 0(6^^^'^). The energy 
level coefficient Ei then approaches infinity as 5 — )■ 0, and thus for small 5 one has that Ei >> 
E^°^. This means that the presence of the SH potential has strengthened the Coulomb interaction 
of the particles, effectively increasing repulsion between the particles. 

NUMERICAL SIMULATIONS 

Simulation of eigenfunctions and comparison with free space. The results above may be 
illustrated using numerical simulations. For Fig. 3 we compute the effective field ip^-^-^ for the 
Dirichlet eigenfunctions in ball B^ where in the ball have in B2 the SH potential we use p = 0.01, 
L = 2tc, E = a and inside the cloak we have a potential Qo represented in the form 

Qo{.r) = riXr<si + r2Xsi<r<s2- (72) 
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using parameters si = 0.6, S2 = 0.8, r2 = —50. Using Matlab, we found the value ri = 12.9016 
corresponding to the SH potential. In the figures we will visualize the effective field 

^ -^^(x) = < _ (73) 

/3m(0)$(x), forxG^i, 



defined in formulas ( [42| ) and ( [49| ), see also the formula ( [46| ). 

Let -Ra = 3 and A^mpty be the event "the particle in the empty ball is in the layer {Ra < 
\x\ < L}" and Ash be the event "the particle in the ball with the SH potential is in the layer 
{Ra < \x\ < L}". Then 

n^empty) = 0.5021, F{Ash) = 0.1355. 

Additionally, let Bempty be the event "the particle in the empty ball is in the layer {2 < \x\ < L}" 
and BsH the event "the particle in the ball with the SH potential is in the layer {2 < \x\ < L}", 
i.e., the particle is outside the cloaking device. We obtain 



empty ) 



0.7196, F{Bsh) = 0.1941. 



The conditional probabilities that the particle is in {Ra < \x\ < L}, conditioned on it being in the 
layer {2 < |x| < i?/,}, are given by 

F{A,,r.pty\Bempty) = n^SHlBsH) = 0.6977. 

This shows that when the SH potential is inserted in the ball B^, the particle is observed in the 
region B2 — -Bi with a lower probability but, when it is observed, the observations from it are 
similar to the observations one would have if the ball were empty. 

Numerical simulation of the scattered wave. We compute scattering solutions for the SH 
potential using spherical harmonics n < A^ = 70 and \m\ < n. In Fig. 1 we plot the total 
field corresponding to the plane wave with energy E = 256 which scatters from the SH potential 
supported in the ball B2. The figure shows the real part of the effective field in the ball B^, L = 3 
in the plane z = 0. The SH potential corresponds to the parameter p = 0.01 and the potential Qq 



is represented in the form (72) with si = 0.25, S2 = 0.5, T2 = —10. Using Matlab, we found the 
value Ti = —169.49 corresponding to the SH potential. 

Simulation of different modes of the cloak. We next visualize the above choice of parameters 
using numerical simulations of the scattering problem. We compute the real part of effective field 
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corresponding to the total field when we have an incoming plane wave, where in the ball B2 there 
is a potential consisting of a cloaking potential V^ ^ plus a potential Qq supported in the cloaked 
region. We use energy E = A, and for the cloaking potential we use the parameter p = 0.01 



and the potential Qq is given in the form (72) with paremeters si = 0.6, S2 = 0.8, T2 = —25. 



Using Matlab, we found the value ri = 7.0675 corresponding to SH potential. We computed the 



corresponding effective field, ip^H' defined in formula (73). To show different modes of the cloak, 



we perturbed the parameter ri to the value rf""-'' = 8.2531 corresponding to the cloak mode and 
computed the corresponding effective field ip^ilak- ^^ addition, we perturbed the parameter ti to the 
value value r^^^ = 7.120 corresponding to the resonance mode and computed the corresponding 
effective field V^^e/. All fields are computed using spherical harmonics n < N = 30 and \m\ < n, 
and the real parts of the fields are plotted in Fig. 2 on the positive x-axis {(x, 0, 0) : x G [0, 3]}. 

IMPLEMENTATION OF THE SH POTENTIALS 



Equation (33 ) with isotropic mass rup^s and bulk modulus k/j describes an approximate acoustic 
cloak. The negative values of T2 required in our construction of potential Qp corresponds then to 
a material with a negative bulk modulus; such materials have already been proposed [!37l . There 
are many proposals for acoustic cloaks ^^ |20l; implementing an acoustic cloak and placing 
negative bulk modulus material inside the cloak, one could test the concept of the SH potential in 
the acoustic setting. 

For electromagnetic waves, one can consider a cylindrical electromagnetic cloak BH and insert 
in it material with negative permittivity to implement a structure similar to the SH potential. The 
analysis related to such cylindrical cloaks with suitable chosen parameters to create a SH potential 
for incident TM-polarized waves is similar to the analysis for the 3D cloak considered in this paper 



(with different asymptotics of coefficients (30)). 

Next we consider approximate quantum cloaks and Schrodinger hat potentials. 

To implement an anisotropic Schrodinger hat, one could consider a quantum cloak with an 
anisotropic effective mass and include the potential Qp in the cloaked region. Possible realizations 
of quantum cloaks have been proposed by Zhang et al. ^, using crystal structures in an optical 
lattice at ultra-low temperatures, which make possible a large variation of the effective mass. 
Below, using the Liouville gauge transformation ([37]), we propose a realization of the SH potential 



using solid state models which do not require large variation of the effective mass. 
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The function V^ ^ is a rapidly oscillating radially symmetric potential, and we can use the ho- 
mogenization theory to approximate the SH potential Vp ^ + Qp by a piecewise constant radially 
symmetric function, 

N 

i=i 
where Xj are indicator functions of suitably chosen intervals /• = [a(j), b{j)] and Vj are con- 
stants. Applying homogenization theory one can see that the solutions corresponding to Vsh with 
sufficiently large A^ approximate the solutions corresponding to Vp^^ + Qp when Vj are chosen ap- 
propriately. In fact, it is enough to use potentials for which the constants Vj have only two different 
values, one very negative and one very positive. Thus, we can form an approximate SH potential 
using layers of spherical potential wells of depth —V- and barrier walls of height V+, i.e., 

VsHir) = 5^ VVx/drl) - ^ ^"^7(1^1) 

where xf ^e indicator functions of the suitably chosen intervals if" = [a^(j), &^(j)]. Note that 
the two-scale homogenization theory used for the above approximation does not specify how large 
A^ one needs to use, but just states that the convergence to the correct limit happens as A^ grows. 

Let us choose the physical units so that mo = 1 is the effective mass of the particle which we 
consider and that h = 2^/^. If ^ satisfies 

. ^2 



-V-V + Vs^(r)-E)^(r)=0, 

then by scaling the length variable by i, we see that \l'(r) = ^(r/£) satisfies 

(_ JL V . V + r'^VsHir/i) - r2E)vl>(r) = 0. (74) 

Zmo 

When i is very large, this means that the energy level E is replaced by a new, much smaller energy 
level i~'^E and the depth of the potential wells are replaced by much shallower wells, and the 
height of the barrier walls are replaced by much lower walls. We denote VsH,ei^) = ^~^KsH(r/-^) 



and Eg = i E. When the above potential Vp^s + Qp is defined using energy E in formulas (39) 



and (40), we say that Vsh/ is designed to operate at the energy level Eg. 

To consider models appearing in solid state physics, we have to consider an effective mass de- 
pending on the X variable. So, now consider four materials with isotropic (i.e., spherical) effective 
masses 1711,1712,171^,7714, such that rrii < ??i2 < rn^ < m^ and rrii < itiq < 777,4 and potentials 
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yi,V2,V3,V4, corresponding to the conduction band edge energies. Assume that the maximum of 
Vj is larger than the maximum VsH,e and the minimum of Vj is smaller than the minimum VsH,i- 
We next study spherical layers of these materials, with abrupt interfaces between materials. 

Let us consider a structure consisting of many thin spherical layers of these four materials which 
near |r| = r are mixed according to ratios ii{r) , i2ir) , (^{r) , i^^r) G [0, 1], correspondingly. We 
need these ratios to satisfy the system 

4(r) + £2(0 +4(r)+ 4(0 = 1, (75) 

mi£i{r) + m2i2ir) + m3i3{r) + mdi{r) = mo, 

rrii 7712 "^3 7714 ^^0' 

VMr) + V2i2{r) + VMr) + V^r) = VshA^)- 



Below, we assume that mj, Vj, Ei, and £"^ V± are such that the solution of the system (75 1 satisfies 
4('^), 4(^), 4('^), 4(^) ^ [0, 1]. By means of homogenization theory, one can see that when the 
equations fTSj ) hold, we can approximate the SH potential using a configuration where spherical 
layers of materials (mj, Vj) are combined at the energy level E^, as is shown below. Note that if 
the masses 771 j are equal to mo we need only two materials and approximate the SH potential with 
spherical quantum wells of given depth —i^'^V- and walls of height i^'^V+. In the general case 



(-^V-— V + F,-E,)^(r) = 0, xeL„ i = i{j) (76) 



when the masses are different, we need four different materials to solve equations (75 1. 
Now consider the situation where we have spherical layers 

Lj = Sk(,) - Sr(,_i) C M^ J = 1, 2, ... , 4 J 

with < R{j) < R{j + 1) < 2. Let Tp be the wave function corresponding to the particle in this 
layered structure. Then in the each layer we have Schrodineger equation 

2 mj 

and i{j) E {1,2,3,4} indicates which material is present in the layer Lj. We may choose i(j) = j 
mod 4 and 

i?(j + l)-i?(j) = ^^.0)(i?(j))- 

On the interfaces of the layers, that is, at r = R{j), we impose the BenDaniel-Duke boundary 
conditions P2ll 

7p\r=R(j)+ = V'|r=R(j)-, drllj\r=R{j)+ = dr^/j\r=R{j)-- (77) 

miQ-+i) mi(j) 
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Similar Hamiltonians and the appropriate boundary conditions in various heterostructures are dis- 
cussed extensively in the reference [|38l . When number of layers 4 J >> A^, the solutions of the 
obtained equations approximate the solutions of the equation f74] ), on related mathematical theory, 
cf. [1231133. 

We now propose how the above model could be physically implemented using effective- 
mass theory to approximate electrons in semiconductors. Consider a semiconductor heterostruc- 
ture build up using four semiconductor materials having the same lattice structures, such as 
Alj.Gai_2,As with mixing parameter x having the values Xi, X2, ^3, and X4 where < x^ < 1, and 
assume that the lowest energies in each conduction band (the conduction band edge) in these ma- 
terials corresponds to the wave vector ko = 0. Using Bastard's envelope function approximation 
[|29ll40ll . we consider in the hetero structure the wave function V^(r), corresponding to the energy 
8 being close to a conduction band edge, which in each material can be expanded as 

V^(r) = 5^/„(rK(r,ko) (78) 

n 

where the sum in n is taken over the finite number of energy bands. Here, /n(r) are the slowly 
varying envelope functions and M„(r, ko) are the periodic Bloch functions in the material corre- 
sponding to wave vector kg. We also assume that the Bloch functions are the same in all four 
materials. Next we consider the single band analysis of electrons near the lowest conduction band 



energy and study the wave function ^(r) omitting in the sum (78 1 all other values of n except the 
value no corresponding to the lowest energy in the conduction band. This means that we use the 
approximation 

^(r) = /„„(r)w„„(r,ko). (79) 

For the GaAs-Ga(Al)As heterostructures, with sufficiently thick GaAs layers, the above ap- 
proximation fTQ] ), with the Schrodinger equations 

(-yV-^V + \/,-^)/„„(r) = 0, rGL,, i = iU). (80) 



and the BenDaniel-Duke interface conditions (77) for the single band envelope function fnai^)-, 



have been proposed in [1291 Sec. 3.II.2.3] and [[4T][. Similar models have also been proposed for 
other 3D heterostructures in [[38l[42l[ . 

Now suppose that the spherically layered semiconductor structure described above is located in 
the ball B2, surround it with a similar heterostructure having effective mass niQ and the conduction 
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band edge energy Vq, being normalized to have the value Vb = 0. We note that above we could 
have assumed that, e.g., 1112 = rriQ and V2 = in which case the surrounding material may 
be homogeneous semiconductor material. To consider the scattering of electrons traveling in the 
medium having the energy 8, one can study the Dirichlet-to-Neumann operator A : h ^-^ d^fno \dB2 
for the boundary value problem, 

(-T^-^^ + ^W)^"oW = ^^"oW' forrG52, (81) 

2 m[r) 

/no(r)|9B2 = ^5 



see (55)-(56). Here, m(r) = m^j-j and V{r) = V^q) in spherical layers Lj C -82- 



As the number of layers, J, grows, the envelope function /„(, (r) approaches [j23l |39l the solu- 
tion of the Schrodinger equation (—2^^ ' ^ + ^SH,e — £)f = ^ and, moreover, the Dirichlet-to- 
Neumann map A of the equation ( (8T| ) approaches the Dirichlet-to-Neumann operator correspond- 
ing to the equation 

;?2 1 

V ■ — V/'^°'"(r) = Sf°^{Y), for r e B^ 

corresponding to the homogeneous background. Thus, the scattering of electrons caused by the 
heterostructure in B2 is very small but the wave function may be very large inside the ball B^. We 



emphasize that Vsh depends on the energy level Ei, see ( |39| ) and ( [40| ), and thus the above analysis 
applies only for electrons whose energy E is close to Ei. 

We note that the above theoretical model can be considered as a (much) more complicated 
structure than the spherical semiconductor layer construction previously used to implement quan- 
tum dots [|43l and a related construction of cylindrical semiconductor layers [|44l . 

Let us also discuss the distribution of the energies of electrons in the heterostructure with con- 
duction band edge energy E^,. The density function of the energies of the electrons in the conduc- 
tion band is the product 

n{E) = D{E)f{E) 

of the density of states D{E) and the Fermi-Dirac distribution f{E) = C'(e(^"^)/(''sT) _^ j^^-i^ 
where T is the temperature, n < Ecis the Fermi energy, ks is Boltzmann's constant, cf. [|45l Sec. 
1.4] and [ l46l Sec. 8], and C is the normalization constant. The density of states D{E) is usually 
approximated by C(£^ — £^c)^^^5 ^ > EcnearEc- When the Fermi-Dirac distribution is approximated 
by Maxwell-Boltzmann distribution fmbiE) = Ce^^^^'^^^^'^^'^K one sees that the density of the 
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energy 8 — 8ch approximately distributed according to the Gamma distribution with the shape 
parameter 3/2 and the scale parameter ksT. Then the energies have the expectation 

3 

and the variance |(A;bT)^. Thus, the energies near the conduction band edge Sc are crucial in 
the modeling of semiconductors, and at low temperatures T the density of electrons, n(£), is 
concentrated near the average £av 

By the above, when the temperature T is low enough, then the energies of the electrons in the 
conduction band have a distribution concentrated near the energy level £av Thus by perturbing 
a homogeneous semiconductor material by including in it one or several SH potentials Vsh,i, 
designed to operate at the energy level Ei = £av, one can create a device where most of the 
electrons would behave in the above analysis: The non-normalized wave functions of electrons 
would not be perturbed outside the supports of the SH potentials but the amplitude of the wave 
functions are strongly amplified inside the support of the SH potentials. 

In summary: Assume that we have layers of semiconducting materials where the electrons with 
energy near the edge of the conduction band can be modeled by a Schrodinger equation ( [8T| ) and 
that in these semiconducting materials the effective masses and the potentials have slightly dif- 
ferent values. Then, appropriate choices of layer thickness and a sufficiently large hetero structure 
lead to the envelope functions of the wave functions satisfing the Schrodinger equation for a SH 
potential. 
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